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Abstract 

QJ We examine the time-dependence of the position-momentum and position- velocity uncertainties 

j^ ■ in the monolayer gapped graphene. The effect of the energy gap to the uncertainties is shown 

^ , to appear via the Compton-like wavelength Ac- The uncertainties in the graphene are mainly 

contributed by two phenomena, spreading and zitterbewegung. While the former determines the 
^ ■ uncertainties in the long-range of time, the latter gives high oscillation to the uncertainties in the 

f^ . short-range of time. The uncertainties in the graphene are compared with the corresponding values 

^: 

lO \ for the usual free Hamiltonian Hjree = (Pi +P2)/2-^^- It is shown that the uncertainties can be 

f~^ ■ under control within the quantum mechanical laws if one can choose the gap parameter Ac freely. 



X: 



I. INTRODUCTION 



)hene(l|, there are a lot of 

n 

activities for researching into the various properties of graphene|2|. This is mainly due to 



After success for fabricating the monolayer or few layer grap 



the fact that the low-energy electrons in graphene have unusual electronic properties. 

Long ago it was predicted by Wallace |3| that the electron located near the hexagonal ver- 
tices of the Brillouin zone exhibits a linear dispersion relation and 40 years later SemenoffJ4| 
showed that the low-energy dynamics of the corresponding electron is governed by massless 
Dirac equation even in the non-relativistic regime. Thus, the fabrication of the monolayer 
graphene opens a possibility to test various predictions of quantum electrodynamics (QED) 
by making use of condensed matter experiment. However, this does not mean that all phe- 
nomena QED predicted can be realized in the graphene-based experiment because the light 
velocity c in QED should be replaced by the Fermi velocity vp ~ c/300. It results in the 
large fine-structure constant a ~ 2. This fact implies that only non-perturbative characters 
of the planar QED can be realized in the graphene experiment. Recently, there have been 
many efforts on this connection |5|. 

Among many phenomena arising in the planar QED most interesting issue, at least for us, 

P 

is the spin-1/2 Aharonov-Bohm (AB)[6] or Aharonov-Bohm-Coulomb (ABC) problem, which 
was extensively discussed about two decades ago[7| because the same problem appeared in 
the context of anyonic and cosmic string theories 8(] . The most important issue in this 
problem is how to treat the 5-like singular potential generated by an interaction between 
particle's spin and thin magnetic flux tube. Recently, similar AB and related problems were 
discussed theoretically |9| and experimentally 10| in the branch of graphene physics. Another 
closely related issue in the graphene is Coulomb impurity problem[ll|. The interesting fact 
in this case is that depending on the charge of impurity there are two regions, subcritical and 
supercritical, in which the effects of impurity are completely different. Similar phenomenon 
in QED was discussed long ago in Ref.|l2[. 

Other unobserved interestin g ph enomena QED predicts are Klein paradox and zitter- 



bewegung. The Klein paradox ll 31 -counterintuitive barrier penetration in the relativistic 

n n 

setting- was re-examined in Ref . \\M, ■ Authors in Ref . [Ij] argued that the Klein paradox can 
be realized using electrostatic barriers in single- and bi-layer graphene. Few years later it was 
reported that the Klein tunneling was observed by measuring the quantum conductance os- 



cillation and phase shift pattern in the extremely narrow graphene 15|. The zitterbewegung 



(ZB)[l6!]-the trembhng motion arising due to the interference between positive and nega 



tive energy states- was also investigated recently in the graphene without |l7l] and with[l8| 
external magnetic field. The effect of zitterbewegung for other models is also discussed 
recently. [19|] 

Besides a connection between graphene and QED much attractive attention is paid to the 
graphene as a new material for future technology. Most important application of graphene, at 
least for us, is a possibility for realization of quantum computer. Recently, many techniques 
are used independently or cooperatively to realize the quantum computer. The typical 
techniques are optical ones, ion traps, NMR, quantum dots, and superconductors. Current 



status for the realization is summarized in detail in Ref. 20|. Also, the graphene-based 
quantum computer is explored in Ref. [211] . 

In this paper we will examine the position-momentum and position- velocity uncertainties 
of the low-energy electrons in the monolayer gapped graphene when the initial wave packet is 
chosen as a general Gaussian wave packet. Since Gaussian wave packet, in general, contains 
both positive-energy and negative-energy spectra, the expectation values of the physical 
quantities should be contributed by spreading and zitterbewegung phenomena. Thus, it is 
of interest to examine the effect of the gap parameter in the expectation values of various 
quantities and uncertainties. We will show in this paper that the position-momentum and 
position-velocity uncertainties can be under control within the quantum mechanical laws if 
the gap parameter can be chosen freely. 

Although this controllability of the uncertainties is an interesting fact from the aspect of 
purely theoretical ground, it is also important from the aspect of quantum computer. Quan- 
tum computer 20| is a machine, which performs quantum computational processes by making 
use of the quantum mechanical laws. So far many quantum information processes are devel- 
oped such as quantum teleportation 22| , factoring algorithm 23| , and search algorithm |24|. 
All quantum information processes consist of three stages: preparation of initial states at 
initial stage, time evolution of quantum states via various unitary gates at intermediate 
stage, and quantum measurements at final stages. If uncertainties, therefore, are large at 
the final stage, the quantum measurement can generate fatal errors in the computing pro- 
cesses. For this reason it is important to reduce the uncertainties as much as possible at the 
final stages. 



This paper is organized as follows. In section II we examine the position-momentum 
uncertainties in the gapped graphene. It is shown that the uncertainties are contributed 
by the spreading and ZB effects of the given wave packet. The uncertainties in the gapped 
graphene are compared with the corresponding quantities of the 2d free Hamiltonian system. 
In section III we discuss on the posit ion- velocity uncertainties in the gapped graphene. 
Unlike the position uncertainties the velocity uncertainties are shown to be contributed by 
sorely ZB effect of the wave packet. This fact implies that the t — )■ oo limit of the velocity 
uncertainties coincides with the Fermi velocity vp regardless of the choice of the packet. In 
section IV a brief conclusion is given. 

II. POSITION-MOMENTUM UNCERTAINTY 

In this section we examine the position-momentum uncertainty in the gapped graphene. 
The appropriate Hamiltonian for the low-energy electron near the Dirac point is given by 

Hm = vf \ (1) 

where vp ~ c/300 is a Fermi velocity and M is a gap parameter generated via some dynam- 
ical and technical reasons. Theoretically, most popular mechanism which generates the gap 



is a chiral symmetry breaking[25]. This mechanism is similar to dynamical breaking 26| . 
which was studied deeply in gauge theories. The bandgap can be generated by breaking 
the sublattice symmetry. This case was experimentally realized by choosing the substrate 



appropriately 27| . In addition, the gap is also generated in graphene nanoribbon 28| . Both 
cases are taken into account in Hamiltonian ([1]). Although monolayer graphene itself does 
not have a gap, the bandgap is naturally generated in bilayer graphene(29|. However, we 
cannot use the Hamiltonian ([T]) to explore the effect of gap in the bilayer graphene due to 
non-trivial structure of the gap in the bilayer system. From the terminology of relativistic 
field theories this gap parameter M is a mass term of the Dirac fermion. 

The position operator x{t) in the Heisenberg picture can be expressed by 2 x 2 matrix 
from x{t) = exp{iHMt/h)x{0) exp{—iHMt/h). Explicit calculation shows 

.w^*,o).if' ^' /''^';-<'"i, (2) 

ai{p)-ta2{p) -S(p) 



where 

0-2 (P) 



h 



p2 + (Mfi.)2 

h 



P2 sin^ 9m + 



{Mvf)pi 



sin Om cos ^ 



A/ cub UM 



v/p2 + (Mt;^)2 ■■ 
[^A/p2 ^ gjn ^^^ pog ^^^ |p2 ^ (Mt;^)2}] 

P1P2 (sin Om cos ^a/ - Om) + (Mt;^)^/^^ + (Mt;^)2 sin^ ^ 



(3) 



[p2 + (Mt;^)2]3/2 

and ^Af = {vFt/K)^Jp'^ + (Mt>i?)2. Each operator in Eq.([n]) consists of two kinds, one of 
which is responsible for ZB phenomena and the other is responsible for the spreading of 
wave packet. 

In order to examine the uncertainty relations we should introduce a wave packet. In this 
paper we introduce a usual two-dimensional Gaussian wave packet 



(x,i/:0)) 



d 



27^^/7^ 



Sk exp 






Ak-r 



(4) 



where real parameters a and h satisfy a2 + 52 _ ]^_ j^ jg easy to show that \ip{x^y : 0)) can 
be decomposed as 

|V^(x, y : 0)) = \r{x, y : 0)) + \r{x, y : 0)), (5) 

where \il)'^{x^y : 0)) and \il)'^{x^y : 0)) are the positive-energy and negative-energy com- 
ponents of \ip{x,y : 0)), respectively. Using Hamiltonian Hm it is easy to derive these 
components and the explicit expressions are given by 

d 



\r{x,y:0)) 



a 



d fcexp 



r/2 ^2 



Ak-r 



(6) 



X- 



k+ + bj^/W+K^ - K') / v/fc2 + A-2 + A-i 

k+^k^ + X-^ I k^ 



\r{x,y:0)) 



d 



d^k exp 



r/2 r/2 



A^T^/7T 

ak+-h{^WT)^ + \-^) I ^/}^T\, 



Ak-r 



A; 



X 



^fc2 + A-2 



-A;. 



In Eq. ([6]) fc-t = k^^iky and Ac = h/lMvp). The parameter Ac is a familiar quantity. In fact, 
this is a Compton wavelength if the Fermi velocity vp is replaced with light velocity c. In 
this paper we will call Ac as Compton wavelength. Thus, the intensity for the positive-energy 



and negative-energy components are 

P+ = (^P(x,y:0)|^P(x,2/:0)) 
P_ = (^"(x,y:0)|V"(x,y:0)) 



2 
1 



AP 



AP 



(7) 



P^ 



where 



AP 



2n 



[d^kexp [-d'iK - «)' - d% - /3)2] ^c'(«'_g^ 
J Wk'^ + , 



2abkj 



(8) 



If, therefore, a = with a = b = 1/a/2, we get P+ = P_ = 1/2. In this case the 
expectation values of various operators are summarized in Appendix A. For arbitrary a and 
/3, however, P± should be computed numerically. Since |'?/'(x, y : 0)) has both positive-energy 
and negative-energy components, the expectation value of various physical quantities should 
exhibit the trembling behavior due to the interference of these components as discussed in 
Ref.[16|-|l9|. 

Using Eq.([2]) and Eq.(jl]) it is straightforward to show 

(a;)(t) = (^(x,2/:0)|x(t)|V^(x,i/:0)) 



TT 



fd'kexp [-d\K - a)' - d\ky - I5f] [Xs + Xzb) 

(9) 



where 
Xs = 

Xzb 



{vFt) 

" fc2 + Ar2 



[{a^ -h'^)\-^K + 2ahk'^^ 



k„ sm U : sm t/ cos y 



'k^ + A; 



+ 



2ab 



(fc2 + A-2)3/2 



(10) 

sin^cos^(A;2 + A;2) 



and, 6 = {vFt)^/k'^ + A"^. As remarked before Xs and Xzb are responsible for the spread- 
ing and trembling motion in the time evolution of the packet, respectively. It is worthwhile 
noting that the fc- integration in Eq.Q can be performed explicitly by making use of the 
binomial expansion. Finally, then, (x) (t) is represented in terms of the Hermite polynomi- 
als. Instead of integral representation, however, (x) (t) has triple summations. The explicit 
expressions in terms of the Hermite polynomials for various expectation values derived in 
this paper are summarized in Appendix B. 

Similar calculation procedure derives {y){t) as 

{y){t) = {^{x,y: 0)\y{t)\^ix,y : 0)) = ^jd^kexp [-d\k, - a)^ - d\ky - pf] {Ys + Yzb) 

(11) 



where 

Ys = -^^, [(«' - b')X-% + 2ahKky\ (12) 

Of course, Ys and Y^^ represent the spreading and ZB motion of the wave packet in y- 
direction. 

In order to confirm the vahdity of our calculation we consider the case of zero gap (A^^ — >■ 
0), which was considered in Ref. [17|. For simplicity, we choose a = 0, a = 1, and 6 = 0. 
Then, Y^ = and Yzb = — sin^^fca,/k^, which makes {y){t) = due to A;a.-integration. In 
this case we also get Xs = and Xzb = sin^ 6'A;y/k^. Using f^^ d9 sin 9 e"'^^'^^ = 27r/i(a), 
where Iu{z) is a modified Bessel function, one can show directly 

{x){t) = ^ (l - e-^''') - de-^''' j^ dqe-'^' cos {^-^q^ h{2/3dq), (13) 



which exactly coincides with the second reference of Ref. 17| . 

Before we explore the uncertainty properties it is interesting to examine the limiting 
behaviors of (x)(t) and {y){t). In the t — )■ limit some combinations of the spreading and 
trembling motion become dominant and the limiting behaviors reduce to 



lim(x) (t) = 2ab{vFt) + O {{vpt)'^) (14) 

lim(y)(t) = {vptf [-{a^ - b^)a + 2abX-^] + O ((M^) • 



It is interesting to note that the t — )■ limiting behaviors of (x) (t) and (y) (t) are completely 
different because their orders oivpt are different from each other. Furthermore, the dominant 
terms of {x){t) come from the off-diagonal components of x{t) while those of {y){t) are 
contributed from all components of y{t). In the t — > oo limit the dominant terms in {x){t) 
and (y) (t) are contributed from spreading terms and their expressions are 

lim (x)(t) = ii!^ Ua^ _ 62)A;Vi,o + 2abJ2,o\ (15) 

lim(y)(t) = ^^^^ [{a' - 6^)A;Vo.i + 2abJ,,,] , 



where 



/k"^k^ 
d'kexp [-d'ik, - a)' - d\ky - /3)2] ^,^^l_r (16) 



In order to examine the position uncertainty Ax{t) we should derive x'^{t), which reduces 
to 

x^{t) = \x\0) + t\p) + ajip) + aiip)] 1 + {x{0),x{t) - x(0)} , (17) 

where {A,B} = AB + BA. Since it is straightforward to show {ip{x,y : 
0)\\ x(0), Z{p) [ \il>{x,y : 0)) = with Z = S, ai, or (3"2, one can show directly 

(18) 



{x'm = y + 7 / d'kexp [-d\K - af - d\ky - Pf] [Xs + Xzb 



where 



X, 



{Vptf 



kl 



fc2 + Ar2 



X.. = sin^^ ^' + ^" 



(fc2 + A- 



2^2■ 



Similar calculation shows 



{y^m = y + 7 / d'ke^v [-d\K.. - af - d\ky - ^Y] (Ys + Yzb 



(19) 



(20) 



where 1^ and Yzb are obtained from Xs and Xzb by interchanging k^ with ky. 

For the case of zero gap (A^^ — ?■ o) with a = 0, a = 1, and 6 = one can show 
straightforwardly 



(.,„4.^,._. 



/32d2 



(21) 



+d'e-^'''' 



d^ 



dq 



Q 



, 2vFt 
1 — cos I q 



{y'm = - + ivFtr 



2^2 



e '^^'^^/^ ( sin ^ + cos ^ ^ 



qWPdq) 
1 



-M'2/3dq) 



d 



-/32d2 



dq 
? 



, 2t;^t 
1 — cos I — —q 
d 



2 y 2/32rf2 
h{2Pdq). 



2(3d 

1 - e-^'"' 



Eq. flT3l) and Eq. fl2T]) can be used to compute the uncertainties Ax and Ay for the case of 
zero gap. 

While in the t — )■ limit {x'^){t) and {y'^){t) exhibit similar behavior as 

lim(x2)(t) = lim(y^)(t) = ^ + {v,tr + O {{v,tr) , (22) 

and the t — )► oo limits of (a;^)(t) and {y'^)(t) reduce to 

(p (p (p (p 

lim (x2)(t) = — + —{vFtfJ2fi lim (i/2)(t) = + _(^;^t)2jo,2. (23) 

Since it is easy to show Ap^ = Apy = h/ \/2d, we plot the time-dependence of the 
dimensionless quantity AxAp^/h in Fig. 1. In the figure we choose a = 0.9, d = 8(nm), a = 



8 



^ 0.503 





Fig. la 



Fig, 1b 



t(fs) 





Fig. 1c 



Fig. Id 



}. (1/nm) 



FIG. 1: (Color online) The time-dependence of AxApx/h for A~^ = 6{l/nm) (a), A^^ = 2{l/nm) 
(b), and A^^ = 0.14(l/nm) (c). The black solid line for each figure is a corresponding value 
{AxApx/h)free for the usual two-dimensional free Hamiltonian Hfree- As Fig. (a), (b), and (c) 
show, the uncertainty AxApx in graphene is larger (or smaller) than {AxApx)free in the entire 
range of time when A~^ > /X2 (or A~^ < /ii). When ^i < A^^ < H2: AxApx is larger and smaller 
than {AxApx)free at t — 7- and t — t- oo limits, respectively, (d) The critical value /i2 increases with 
decreasing a, and eventually goes to oo at a = 0. 



0.04(l/nm), and /3 = 1.2(l/nm). We also choose the inverse of the Compton wave length 
as Q{l/nm) (Fig. la), 2{l/n'm) (Fig. lb), and 0.14(l/nm) (Fig. 1 c). The black solid line 
in (a), (b), and (c) is {AxApj,/K)free = a/(1/2)2 + {\cVFt/2d?Y, which is a corresponding 
value for the usual non-relativistic free Hamiltonian Hfree = {Pi + Pq)/'^^- The unit of the 
time-axis is femto-second. 



As Fig. 1 represents, the uncertainty AxAp^ has several distinct properties. First, it is 
contributed from both spreading and ZB motion of the wave packet. The spreading motion 
is dominated in the large scale of time. With increasing the inverse Compton wavelength 
the overall increasing rate of AxAp^,' arising due to the spreading of the packet decreases 
drastically. This can be understood from the analogy of the relativistic field theories, that is, 
with increasing M the relativistic theories approach to the non-relativistic Galilean theories, 
where the uncertainty is minimized. In the small scale of time Ax^px oscillates rapidly due 
to the ZB effect. The amplitude of the oscillation increases with decreasing A~^. This is 
mainly due to the fact the the ZB effect is dominated when the energy gap /\E between 
positive and negative energy spectra decreases. However, the frequency increases rapidly 
with increasing A~^ because of the famous formula u = AE/h. When A~^ is larger than 
a critical value /i2, AxApx becomes larger than [Ax Ap^) free as Fig. la indicated. When, 
however, A~^ is smaller than a different critical value /ii, it is smaller than {AxApx)free 
as Fig. Ic shows. In the intermediate range of A~^ AxApx is larger and smaller than 
{AxApx)free in t — )■ and t — )■ oo limits, respectively as Fig. lb shows. Using Eq.( IT^ . 
( IT5|) and several other limiting values, one can derive the critical values /ii explicitly, and 
/i2 implicitly as 

1, (24) 



/^l = /^ ,o.. . ^^^^ 7(Ac ^) 



v/2t/2(i_4a262) 
where 



liK') 



'c 



TT 



d f , o ,<,x , _i , „ , , ^2 



Ac =M2 

(25) 



J2.0 - - {(a' - &')A;Vi,o + 2abJ2flY 

71 

The A~^-dependence of 7(A~^) is plotted in Fig. Id when a = 0.9, d = 8{nm), a = 
1.2/n{l/nm), and /3 = 1.2{l/nm) for various n. As this figure indicates, the critical value 
^2 increases with increasing n, and eventually fi2 = oo when a = 0. 

The dimensionless uncertainty AyApy/h is plotted in Fig. 2 when a = 0.9, d = 8(nm), 
a = 1.2{l/nm) and (3 = 0.04(l/nr7i). We also choose A~"^ as 8{l/nm) (Fig. 2a), 2{l/nm) 
(Fig. 2b), and 0.08(l/nm) (Fig. 2c). We plot {AyApy/h)free together for comparison. As 
Fig. 2 shows, AyApy exhibits a similar behavior with AxAp^. However, the critical values 
fj,i and fj,2 are changed into z/i and z/2, which reduce to 

= 1, (26) 
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FIG. 2: (Color online) The time-dependence of AyApy/h for A~^ = 8(l/n?7i) (a), A~^ = 2{l/nm) 
(b), and A~^ = 0.04(l/?7-m) (c). The black solid line for each figure is a corresponding value 
{AyApy/h)free- As Fig. (a), (b), and (c) show, the uncertainty AyApy in graphene exhibits a 
similar behavior to AxAp^- However, the critical values fii and //2 are changed into i^i and 1^2- (d) 
The critical value 1^2 increases with decreasing /3, and eventually goes to 00 at /? = 0. 



where 



6{\ 



2M 



-l^ 



2X-/d 



IT 



\2--{{a'-b')\-'Jo,^ + 2abJ^,,Y . (27) 

TT 

The A~^-dependence of 5(A^^) is plotted in Fig. 2d when a = 0.9, d = 8{nm), a = 
1.2{l/nm), and f3 = 1.2/n{l/nm) for various n. As this figure indicates, the critical value 
1/2 increases with increasing n, and eventually goes to 00 when /3 = 0. 



11 



III. POSITION- VELOCITY UNCERTAINTY 



In this section we discuss on the position- velocity uncertainties 30|, which is completely 
different from position-momentum uncertainties because of p 7^ Mv. The velocity operator 

(0) = dHuldpi. This 



Vx{t) is defined as exp ( i^M^^) "^^lO) exp ( — ii^TV/V^)? where v,^ 
operator is easily constructed from x{t) by making use of Ehrenfest 
(i//i) exp iiHMt/h) [Hm, x{0)]exp i—iHMt/h) = Vx{t). Then, the final expression oi Vx{t) 



31l | theorem dx{t)/dt 



IS 



Vx(t) 



U{p) Ui{p)+iu2{p) 

Ui{p)-iu2{p) -U{p) 



(28) 



where 



U{p) = vf 
ui{p) = Vf 
U2{p) = Vf 



, , , ===■ sm du cos Gm + ^ , .,, — 7^ 



sin^ 6'j\/ 



(29) 



COS^ ^M + ^ 



pi 



(Mf^)^ 



p2 -I- {MvfY 

2piP2 _,__2 



sin^^ 



M 



p2 + (Mt;^)- 



sin^ e 



M 



2{Mvf) 



v/p2 + {MvFy 



sin^M cos^ 



M 



Unlike the position operators x{t) and y{t) the velocity operator Vx{t) does not have the 
spreading term. This is due to the fact that the spreading term in the position operators is 
linear in time. Another remarkable property of Vx{t) is that vl{t) is simply Vp times identity 
operator 1. Combining these two properties one can easily conjecture limt^oo ^Vx = vf 
regardless of the choice of the wave packet because the ZB term in Vx{t) has infinitely high 
frequency in this limit, and therefore, is canceled out in the time average. 

The expectation value {vx)(t) and {v1){t) with a wave packet (jlj) can be straightforwardly 
computed by making use of Eq. fl28l) . As expected the resulting Avx{t) has only trembling 
motion and approaches to fi? at t — )■ 00 limit. The dimensionless position- velocity uncer- 
tainty AxAvx/dvF is plotted in Fig. 3 for A^-*^ = 0.09(l/r2m) (Fig. 3a), A^-*^ = 0.14(l/nm) 
(Fig. 3b), and A""*^ = 0.5(l/n?7i) (Fig. 3c) when a = 0.9, d = 8{nm), a = O.OA{l/nm) and 
P = 1.2{l/nm). The x-axis is time axis with femto-second unit. The black dotted line is 
a corresponding value {AxAvx)free/dvF, where {AxAvx)free = \/^lvp/A + Xf.Vpt'^/Ad'^ is a 
position- velocity uncertainty for Hfree- The overall increasing behavior of AxAf^: is solely 
due to Ax because Avx does not have its own spreading term. As Fig. 3 shows, AxAwa. 
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Fig. 3a 



t(fs) 




Fig. 3b 




Fig. 3c 



FIG. 3: (Color online) The time-dependence of /S.x/^Vx/dvF for A^-"^ = 0.09(l/n?7i) (a), A~^ = 
0.14(l/nm) (b), and A~^ = 0.5(l/n?7i) (c). The black dotted line for each figure is a corresponding 
value {IS.xIS.Vx / dvp) free- As Fig. (a), (b), and (c) show, the uncertainty ISx/Svx in graphene is 
larger (or smaller) than {ISxISvx)free depending on the gap parameter A^^. One can show explicitly 

that Wuit^QlSxISVx < (AxAVx) free if ^c^ < /"I ^l^d liuit^oo AxAVx > (AxAVx) free if ^c^ > /"2*) 

where ^2* is defined as 7(A^^ = fi2*) = l/(2(/X2*(i)^. 

can be smaller or larger than {AxAvx) free depending on the gap parameter Ac. In order to 
compare AxAti^. with {AxAvx)free more accurately we compute its limiting values at t — )■ 
and t — 7- 00. Then, it is easy to show limj^^o AxAf^. < (AxAvx) free if ^c^ < /^i; where /ii is 
defined at Eg. (1211) . and limt_i.oo AxAt>^ > {AxAvx)free if ^c^ > A^2*, where /i2* is defined as 
7(A~^ = fi2*) = l/(2(/^2*'^)^- The critical values /ii, /i2, and ^2* are given at Table I when 
d = Slnm), a = 1.2/n{l/nm), /3 = 1.2{l/nm), and a = 0.9 or 0.7. The reason for choice of a 
is that while the diagonal components of the various operators contribute dominantly to the 
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uncertainty relations at a = 0.9 ~ 1, the off-diagonal components become more important at 
a = 0.7 ~ \l \f^. As expected from Fig. Id, 1x2 increases with increasing n, and eventually 
goes to 00 at q; = 0. Another critical value /i2* also exhibits an increasing behavior with 
increasing n, but its increasing rate is very small compared to /i2 and converges to 0.332 at 
n — > 00 limit. 

Table L Critical values for AxAp^ and lS.xlS.Vx when d = S^nm), a = 1.2/n{l/nm) and 

/3 = 1.2(l/nm). 





a 


n = 10 


n = 20 


n = 30 


n = 40 


n = 50 


n = 00 


/ii(l/(nm) 


0.9 


0.143 


0.143 


0.143 


0.143 


0.143 


0.143 




0.7 


4.42 


4.42 


4.42 


4.42 


4.42 


4.42 


1^2(1 /{nm) 


0.9 


1.03 


2.24 


3.47 


4.69 


5.90 


00 




0.7 


0.90 


1.79 


2.68 


3.58 


4.47 


00 


Id2*{l/{nm) 


0.9 


0.257 


0.303 


0.318 


0.324 


0.327 


0.332 




0.7 


0.256 


0.302 


0.317 


0.323 


0.326 


0.332 



Following similar calculation procedure one can plot the time- dependence of the dimen- 
sionless quantity AyAvy/^dvp)- Although the time-dependence of the uncertainties is not 
plotted in this paper, AyAvy exhibits a similar behavior with AxAf^. However, the critical 
values /xi and fj,2* are changed into ui and z/2*, whose explicit values are given at Table 11. 

Table IL Critical values for AyApy and AyAvy when d = 8(nm), a = 1.2(l/nm) and 

/3 = 1.2/n{l/nm). 





a 


n= 10 


n = 20 


n = 30 


n = 40 


n = 50 


n = 00 


i'i{l /{nm) 


0.9 


0.088 


0.088 


0.088 


0.088 


0.088 


0.088 




0.7 


0.088 


0.088 


0.088 


0.088 


0.088 


0.088 


V2{l/{nm) 


0.9 


2.23 


3.36 


4.48 


5.60 


6.73 


00 




0.7 


1.22 


2.05 


2.88 


3.73 


4.59 


00 


U2*{1 /{nm) 


0.9 


0.309 


0.326 


0.329 


0.330 


0.331 


0.332 




0.7 


0.319 


0.328 


0.330 


0.331 


0.331 


0.332 
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IV. CONCLUDING REMARKS 

In this paper we have examined the position-momentum and position-velocity uncer- 
tainties for the monolayer gapped graphene. We have shown that the uncertainties are 
contributed by the spreading effect of the wave packet in the long-range of time and the ZB 
in the short-range of time. By choosing the gap parameter Ac appropriately one can control 
the uncertainties within the quantum mechanical law. 

The uncertainties can be tested experimentally because all figures in this paper show a 
significant difference between free and graphene cases. The uncertainties in the graphene 
might be measured via the following one-slit experiment (see Figure 4). In this paper we 
will discuss on Ax only because other quantities can be measured similarly. The slit width d 
should be order of Angstroms to ensure the occurrence of diffraction in the slit. The distance 
L should be order of nanometers because the effect of the zitterbewegung is important within 
initial few femtoseconds. The electrons emitted by the emitter would arrive at the detecter 
through the slit. Then, one can make a probability distribution with respect to x, which 
would be a smooth Gaussian form. Measuring the width of the Gaussian distribution, one 
can deduce Ax at t ~ L/vp, where vp is a Fermi velocity. Repeating the same experiment 
with changing L one can measure the time-dependence of Ax. If the prediction we presented 
in this paper is correct. Ax would exhibit an oscillating behavior in the short-range of time 
due to the effect of the zitterbewegung, but globally an increasing behavior in the long-range 
of time due to the spreading effect of the wave packet. 



electron 
emitter 



Graphene Plane 



^ 



detector 



FIG. 4: (Color online) Schematic diagram for measuring the uncertainties. 
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It is interesting to extend this paper to the bilayer graphene. Another interesting issue 
is to examine the uncertainty relations when external magnetic field is applied. We guess 
that the external magnetic field drastically reduce the uncertainties in the graphene. If so, 
the graphene-based quantum computer can be more useful for huge calculations. We would 
like to explore this issue in the near future. 
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Appendix A 

In this appendix we summarize the various expectation values at a = and a = b = l/-\/2, 
where Eq. (JTj) and Eq. ([8]) imply that the initial wave packet has equal intensity of positive- 
energy and negative-energy states. In this simple case the expectation values (x)(t) and 
{y)(t) reduce to 



TT 7 






(A.l) 



71 J K + A„ 



where 9 = {vpt)^y'k'^ + A~^. In the case of zero gap we get {y){t) = 0. Since (a;^)(t) and 
{y'^)(t) are independent of choice of a and b, those are equal to Eq. flTSj) and Eq. fl20|) with 
a = 0. The expectation values for the velocity operators becomes 

{v,){t) = vp-^^ [ d'ke-''''^-''^'y-P^\m'e'i^^ (A2) 

11 J k + A^ 



^ J ^/k^ + A-2 

In the case of zero gap we also get {vy)(t) = 0. Of course, expectation values for the square 
of velocity operators are simply (t>^) = {Vy) = Vp. 
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Appendix B 

In this appendix we summarize the explicit expressions for {x)(t), {y)(t), {x^){t), {y'^){t), 
{vx){t), and {vy)(t) by making use of the binomial expansion and performing the k- 
integration. The integral formula we use is 

^n^-ix-fsf^^ = (2i)-"v^iy„(i/3), (B.l) 

) 

where Hn{z) is the usual Hermite polynomial. 

The expectation values (x)(t) and {y){t) expressed in Eq. ([9]) and Eq. flTTl) reduce to 

(^)(*) - 2*.* ^ t ^^1^ t ( ; ) pg^ (B.2) 

i{a^ - h'^)dXi + 2ah{vFt)X^. 

i{a^ - h')dYi + ahX^y2 







n=0 

where 

Xi = {2n + 3)H2m{iad)H2e-2m+i{iPd) + 2{X~^VFt)H2m+i{iad)H2e-2m{iPd) (B.3) 

-^2 = H2m{iad)H2i-2m+2{i(3d) - {2\^^dfH2m{iad)H2e-2m{i(3d) 

Yi = (2n + 3)H2m+i{i(yd)H2i-2m{if!^d) + 2{\~^VFt)H2m{i0ld)H2l-2m+l{i(^d) 

Y2 = {2n + ?,){2\-;^dfH2m{iad)H2t-2m{iPd) - 2{X;Wt)H2m+i{tad)H2e-2m+i{iM- 

Although the arguments of the Hermite polynomials are pure imaginary, one can show easily 
that (x)(t) and {y)(t) are real by considering the fact that Hn{z) is even (or odd) function 
when n is even (or odd). 

Similarly one can express (x^)(t) and {y'^){t) from Eq. flTSl) and Eq. fl2Ul) as follows: 

{x'm = ^+{vFtr (B.4) 

h (2- +4)! hv I ^^"^-^'^r-^ L 

+''^ (2n + 4)! l.^,j(2A,W-^A.^^r^ 
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where X3 = X2 and 



Y3 = H2rn+2iiad)H2E^2miil3d) - (2 Ac d) H2miiad)H2e-2miil3d). 



(B.5) 



Although we have not derived the integral representations of {vx){t) and {vy){t) explicitly 
in the main text, their derivations are straightforward. Then, the expressions of {vx){t) and 
{vy){t) in terms of the Hermite polynomials are 



{vx)(t) = 2ahvp — 2vp 



00 

n=0 



(2n + 2) 



2n+l " 



n 



\n-i 






i (2A-irf) 



M^2£+2 



(B.6) 



m=0 



m 



i{a' - ¥){2\-^d)Ui + 2ah{\^\Ft)U2 



{vy){t) =vfY^ 

ra=0 



n 



^2A£Vt)^Y 

(2n + 2)! f;^\ £ / {2\-^dr+ 




i{a^ - h'^)Vi + 2ahV2 



where 



Ui = {n+ l)H2rn{iad)H2i-2m+i{i(3d) + {\~^VFt)H2m+iiiad)H2i-2m{i(3d) (B.7) 

U2 = {2\^^dfH2m{iad)H2e^2m{if3d) - H2m{iad)H2i-2m+2{i(id) 

Vi = {2n + 2)H2rn+l{iad)H2t^2rn{il^d) - 2{\~^VFt)H2m{iad)H2t-2rn+l{i(^d) 

V2 = {2n + 2){2\^^d)H2m{ioid)H2e^2m{if3d) —H2m+i{i<^d)H2i-2m+i{i(^d). 
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